Andreev tunnelling in quantum dots: A slave-boson approach 
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We study a strongly interacting quantum dot connected to a normal and to a superconducting 
lead. By means of the slave-boson technique we investigate the low temperature regime and discuss 
electrical transport through the dot. We find that the zero bias anomaly in the current-voltage 
characteristics which is associated to the occurance of the Kondo resonance in the quantum dot, is 
enhanced in the presence of superconductivity, due to resonant Andreev scattering. 
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With the advent of nanotechnology, recent years have 
witnessed an impressive experimental activity, studying 
various properties of small mesoscopic structures. In 
particular, the transport properties of hybrid supercon- 
ducting structures and the associated Andreev scattering 
mechanism have been investigated intensively after ob- 
serving several new phenomena ^ . Many of these phe- 
nomena have been successfully explained in terms of a 
one-particle picture essentially based on the BCS theory 
either via the Bogolubov-de Gennes equations or via qua- 
siclassical Green's function methods, as documented by 
various review articles [Q. 

Besides this, electrical transport through small con- 
fined regions, where electron-electron interactions are 
strong, has also attracted a lot of interest. Such exper- 
imental setups, e.g. quantum dots (QDs), allows one 
to investigate in a controlled way the interplay of the 
electron-electron interaction and disorder. In particu- 
lar it has been pointed out that a QD attached to two 
metallic leads resembles an impurity level in a metal. As 
a consequence, even when the dot level is far from the 
Fermi energy of the leads, transport will occur due to 
the Kondo effect ||]. This is due to the formation of a 
spin singlet between the impurity level and the conduc- 
tion electrons, which gives rise to a quasiparticle peak at 
the Fermi energy in the dot spectral function. This sug- 
gestion has been explored theoretically by several groups 
[^-Q and lead to the prediction that one should observe a 
zero-bias anomaly in the current voltage characteristics. 
Such an anomaly has been indeed observed in different 
QD systems |,§. 

In a recent letter it has been suggested that, if the 
QD is coupled to two different types of leads, i.e., a nor- 
mal and a superconducting lead, resonant Andreev tun- 
nelling yields a stronger zero-bias anomaly with a broader 
temperature region where the effect occurs. In the anal- 
ysis of | pO| an approach based on the equation-of-motion 
method which is mainly valid at high temperature was 
used. In this paper, we extend the analysis of to 



the extreme low temperature regime. To this end we use 
slave-boson mean field theory. This approach has been 
successfully applied to the low temperature properties of 
a Kondo impurity in presence of normal jri|,p^ as well 
as for superconducting conduction electrons |13|. Despite 
its simplicity, this method captures the main physical as- 
pects of the Fermi liquid regime at low temperatures, i.e., 
the formation of a many-body resonance at the Fermi en- 
ergy. For this reason it presents a convenient framework 
in which to study the interplay between Andreev scatter- 
ing and Coulomb interactions. 

We model the N-QD-S system with the Hamiltonian 



H — Hn + Hfi + Hr 



H' 



T,N 



T,S 



(1) 



where Hj^ = Ska ^kcjj ka^N.ka , ^^"3 = 
Eka '^kc^^cs^ka + Ek(^4,kT4,-ki + c-c.) and Hu = 
edd^d^ + Uud^ndi are the Hamiltonians of the normal 
lead, the superconducting lead (A is the superconduct- 
ing gap) and the dot respectively. The single particle 
energy is double degenerate in the spin index a and 
the interaction is included through the on-site repulsion 
U. The position of the dot level can be modulated by an 
external gate voltage. Tunnelling between the leads and 
the dot is described by HT,a = Eko-(^a4,k(T^'^ + ^-C-) 
where a = N, S and Va is the tunnelling amplitude. For 
simplicity we assumed Va independent from k and a. 

In the following, we consider the on-site repulsion U 
is infinite, so processes where the dot level is doubly oc- 
cupied are excluded. This condition allows to apply the 
slave boson technique The dot level is represented 
as dl = f^b, where the fermion and the boson b de- 
scribe the singly occupied and empty dot states. Since 
the dot is either empty or singly occupied, the constraint 
b^b + E<T flu = 1 has to be fulfilled. 

In the mean field approximation, the operator b is re- 
placed by a c- number 60, and the constraint is fulfilled 
only on average. This is achieved by introducing a chem- 
ical potential Aq for the pseudo particles. Notice that one 
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ends up with a non interacting-like problem with renor- 
mahzed parameters, i.e., an energy shift for the dot level 
Ed ^ fid + Ao = id and a multiplicatively renormalized 
tunnelling Vq, —> boVa. 

We discuss the mean field equations and its solu- 
tion first in equilibrium and then generalize to non- 
equilibrium. We start from the impurity part of the free 
energy, which in presence of both normal and supercon- 
ducting leads is given by 

F = -T^ Trln[ie„a" - ?d<T" - blt{ie„)] + Ao6g + - fi, 

(2) 

where e„ is a fermionic Matsubara frequency, ct* are the 
Pauli matrices, and 



f(ie„) 



|V„|^6-^gk,a(ien)<5-'' 



(3) 



with ^k,Q being the Green's function of the lead a. 

By minimizing the free energy with respect to Aq and 
bo we find the equations 



boAo + boT^Tr 



.2+T^Tr[g(i6„)a^^ 
[^(ien)r(ieri) 



= 0, 

= 0, 



(4) 
(5) 



which have to be solved self-consistently. ^(ie„) is the 
pseudo fermion Green function given by ^(ie„) = [iCnO'^ — 
eda"-62f(ie„)]-i . 

Before presenting a numerical solution of the above 
equations, it is useful to get some insight from an approx- 
imate analytical treatment. The first equation, eq.(i|), is 
the constraint. Since the pseudo fermion level is at max- 
imum singly occupied, the renormalized level is above 
the Fermi energy. In the Kondo limit, where the occu- 
pancy is nearly one, it is found < Cd < feo(rN + Ts); 
i.e. Aq ~ |ed| and ed ~ 0. The renormalization of the 
tunnelling amplitude is determined from eq.(^. A triv- 
ial solution bo = always exists. The solutions which 
minimize the free energy, however, are those with bo ^ 0. 
By introducing a flat density of states in the leads and 
the tunnelling rates Ta = ttA'oqIVqP, the elements of the 
matrix r(ie„) are Fn = F22 = — i7i and F12 — T^i ~ 72, 
where 



71 = sign(e„)FN -I- Fs 



p 72 = Ts 



A 



(6) 



Restricting ourselves to zero temperature, we replace the 
Matsubara sum in eq. (|^) by an integral and obtain 



kdl = 4 



w 



de 7i(e + 6§7i)+&ol72|' 



(7) 



where W is a cut-off of order of the band width. We 
simplify the integral by approximating 71 and 72 as 



71 



for e < A 
Fs for e > A 



72 







for 
for 



e < A 
e> A 



(8) 



The result is 



ed 



Fn (A + 62Fn) 



In 



(A + blTs + fegFN)2 ■ 



(9) 



where we neglect a term proportional to Fg , but without 
any logarithmic factor. For small superconducting gap, 
i.e. A much smaller than the Kondo temperature which 
is given by Tk = 6qFn + b^Ts, A is neghgible. One can 
then easily solve eq.(^) for 5q and obtains the result for 
two normal leads with total tunnelling rate Fn -I- Fs : 



6g(FN + Fs) = Tyexp 



ed 



2 Fn + Fs 



(10) 



In the opposite limit, where A is much larger than the 
Kondo temperature, we find 



7r|ed|-(2FsA) HW/A)\ 
2 



(11) 



The results agree qualitatively with what we expect from 
scaling arguments for the Anderson model [Q. In the 
perturbative regime, a logarithmic correction to ed has 
been found. In the case of a large gap scaling due to the 
superconducting electrons stops at energies of the order 
A, giving rise to a finite renormalization of ed, as seen 
in eq.(pl|). A similar shift in the dot level has also been 
reported within the equation-of-motion approach of ref. 
pO| . In case of a small gap, the superconducting lead 
contributes to scaling down to low energies, where one 
enters the strong coupling regime. Presumably, the fixed 
point is still reached for energies of the order of Tk, much 
greater than A, so that the Kondo temperature does not 
depend on A, compare eq . ([lo|) . 

Notice that in presence of normal electrons, we always 
find a non-trivial solution of the mean field equations. 
This is to be contrasted with the case of superconduct- 
ing electrons only, Fn = 0, where for large gap only the 
solution 60 = exists, and there is no Kondo effect 03] . 



2 



1 

0.1 
0.01 
0.001 



0.0001 





1 


1 1 1 1 1 1 1 b 



















\ 




rs = 2rN -- 


N 















0.01 



0.1 



FIG. 1. The mean field parameter bQ, which is a measure 
for the Kondo temperature, as a function of the supercon- 
ducting gap. For both curves the energy of the dot level is 
= -W/3 and the total tunnelling rate Fs + Fn = 0.15W. 
The Kondo temperature estimated from these parameters is 
Tk ^ 0.03W. 

In Fig.|l| we report numerical results for 6q as a func- 
tion of A. Instead of a sharp cut-off in the density of 
states, we used an exponential, exp(— e^/FF^), in our nu- 
merics. The bare d-level is Cd = — FF/3. For Fn + Ts = 
0.15VF the Kondo temperature is approximately given by 
Tk « 0.03VF. As long as the gap remains below this en- 
ergy scale, 6o is almost constant. It drops quickly for 
Tk < A < W, with a region, where log 6g decreases lin- 
early in log A, in agreement with eg. (pi]) . 

In a non-equilibrium situation, when a voltage is ap- 
plied between the two leads, the mean field parameters 
cannot be obtained by minimizing the free energy. How- 
ever one can derive self-consistency equations following 
the arguments given in ref. p^ . These impose the van- 
ishing of so called tadpole diagrams order by order in 
perturbation theory. The equations read 



27r 



-Tr 



g<{e)a' -0 (12) 



Xobo - ibo I ^Tr 

ZTT 



e«(e)F<(e) + g<(e)f-4(e) =0, (13) 



where the lesser Green's function Q^{t,t') — \{(f)^ {t')<j){t)) 
has been introduced, with (j) = {J],j\)- The lesser and 

advanced matrix F is defined in analogy to its equilib- 
rium version m eq. (|). To obtain g< , we use the gen- 
eral relation ^< = Q^Ti^Q^^ where at mean- field level 
t< = 6^f< and 



T<{e) 



E 



Va\ O 



?^k(e)/a(e) 



/a(e)5a,k(e 



(14) 



If the chemical potential of the normal electrode differs by 
eV from that of the superconductor, the matrices have 
elements = /(e-t-eT/^) and /q,22 = l-/(-e-|-eV'^), 
with Vj^ = Vg^ = and /(e) being the Fermi function. 



Note that the superconducting lead does not contribute 
to S<(e) for e < A. 

The solution of the mean- field equations, in presence of 
an external voltage, can then be obtained along the lines 
of the equihbrium case. As long as \eV\ < Tk the solu- 
tion is almost independent of the voltage. For large volt- 
age, |ey| » Tk, we found the Kondo peak pinned to the 
chemical potential at the normal side, i.e. id ^ id — eV, 
with a decreased width. 

Following ref. ]l0| , the Andreev current through an in- 
teracting quantum dot is 



2ie / ^FNTr{a-G«(e) 



S«(6),/n(6)] G^{e)} 



(15) 



where G and S are the Green's function and self-energy 
of the dot electrons. Within the present mean-field ap- 
proach, G = b^Q. Explicitly, the Andreev current is given 
by 

,(^).^ra /'— "-)-/'■ ^-"V hsm (16) 



with 



GNs(e) 



4(fNfs)2 



(6-2 



Here the tunnelling rates Fs,n = ^oi^s.N, 



(17) 



and i — 



6qFs/A). Then one recovers the current formula 



for a non- interacting quantum dot, with renormalized 
parameters which are voltage dependent. We see that 
on resonance, when ~ and e = 0, the small renor- 
malization factor bg drops out and the conductance has 
the same value as in the non-interacting case |18|. The 
peak strength becomes maximal when Fn — Fg with 
GNS,max — 1, twicc the maximum for a N-QD-N system. 

Finally, we want to comment on the reliability of our 
results. The success of slave-boson mean field theory 
stems from the fact that it captures the Fermi-liquid 
regime at low temperature. If the N-QD-S system scales 
to a Fermi liquid at low temperature, Gns as given in 
eq.(|l^) is exact in the low temperature, low voltage limit, 
with the parameters to be determined. Within slave- 
boson mean field theory Fn and Fg renormalize equally, 
although this may no longer be the case when consid- 
ering higher order corrections. For illustration, we es- 
timate the effect of residual quasiparticle interaction in 
the limit A <^ Tk. By assuming an effective quasiparti- 
cle interaction of the form iJint = Un^ni, we find to first 
order in U no corrections to Fn, while, as one could have 
expected, repulsive quasiparticle interaction suppresses 
fg = blVs[l - (i7/7rTK)(A/TK)lnTK/A]. 

We studied Andreev tunnelling through a strongly in- 
teracting quantum dot, focussing on the extreme low 
temperature limit. In agreement with a previous study. 
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we found an enhanced Andreev current at low bias volt- 
age, due to the Kondo effect. The zero-bias conductance 
is maximum with the universal value Gns = 1 when the 
renormalized tunnelling rates Fn and Fs are equal. We 
identified the ratio A/Tk as an important parameter. In 
the case A <C Tk, the Kondo resonance forms as for 
two normal leads. The condition Fn — Fs coincides with 
equal bare tunnelling rates Fn = Fs- In case of large gap, 
quasiparticle interaction suppresses Fs, nevertheless the 
conductance maximum condition may be achieved with 
large bare tunnelling rate Fs > Fn. 
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